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Multi-kink topological terms and charge-binding domain-wall condensation induced 
symmetry-protected topological states: beyond-Chern-Simons/BF field theories 


Zheng-Cheng Gu, 1 ’* Juven C. Wang, 2,1 ’! and Xiao-Gang Wen 2, * 
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Quantum-disordering a discrete-symmetry breaking state by condensing domain-walls can lead 
to a trivial symmetric insulator state. In this work, we show that if we bind a ID representation of 
the symmetry (such as a charge) to the intersection point of several domain walls, condensing such 
modified domain-walls can lead to a non-trivial symmetry-protected topological (SPT) state. This 
result is obtained by showing that the modified domain-wall condensed state has a non-trivial SPT 
invariant - the symmetry-twist dependent partition function. We propose two different kinds of held 
theories that can describe the above mentioned SPT states. The first one is a Ginzburg-Landau- 
type non-linear sigma model theory, but with an additional multi-kink domain-wall topological term. 
Such theory has an anomalous U k ( 1) symmetry but an anomaly-free Z^ symmetry. The second 
one is a gauge theory, which is beyond Abelian Chern-Simons/BF gauge theories. We argue that 
the two Held theories are equivalent at low energies. After coupling to the symmetry twists, both 
theories produce the desired SPT invariant. 
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I. INTRODUCTION 

A. SPT states and their effective field theories 

Recently, it has been realized that many-body ground 
states can be divided into two classes: [1] long-range en¬ 
tangled (LRE) states and short-range entangled (SRE) 
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states. The LRE states can belong to many differ¬ 
ent phases that correspond to topologically ordered 
phases. [2, 3] When there is a global symmetry (described 
by a group G), even SRE states can belong to many 
different phases, and these phases are called symmetry- 
protected topological (SPT) states. [4-9] A large class 
of bosonic SPT states whose boundary has a pure 
“gauge anomaly” [10-12] can be systematically classified 
via group cohomology classes H d+1 (G, R/Z). [13-15] All 
these SPT states can be realized by exactly-soluble lat¬ 
tice non-linear cr-model with the symmetry group G 
as the target space plus a 27r quantized topological 9- 
term. They can also be realized by exactly-soluble lat¬ 
tice Hamiltonians that contain seven-body interactions. 
In addition, bosonic SPT states whose boundary has a 
“gauge gravitational mixed anomaly” can all be realized 
by lattice non-linear cr-model with SO oo x G as the tar¬ 
get space and with a 27r quantized topological 0-term. [20] 
The potentially possible SPT invariants of the first and 
the second classes of SPT states can also be studied di¬ 
rectly via cobordism theory, [16-19] but the cobordism 
theory does not give rise to a realization of the SPT 
states. 

Many of the SPT states protected by discrete group 
symmetry can also be realized by condensing domain 
walls in symmetry breaking states, if we decorate the do¬ 
main walls with lower dimensional SPT states and/or in¬ 
vertible topologically ordered states. [20-23] In this work, 
we will realize some additional SPT states by condensing 
domain walls, such that the intersection point of several 
domain walls carries the quantum number of the unbro¬ 
ken symmetries. More general SPT states protected by 
discrete group symmetry can be obtained by decorating 
the intersection lines (or surfaces) of several domain walls 
with ID (or 2D) SPT states (as indicated by the Kunneth 
formula for the group cohomology [20, 21]). 

In addition to the above systematic constructions of 
all the bosonic SPT phases, people have also developed 
many field theory realizations for some special simple 
SPT states (under the name of bosonic topological insula¬ 
tor (BTI) [23-31, 46]) which lead to some simple physical 
pictures and mechanisms for bosonic SPT states. Due 
to the incompressibility of topological phases, it is suffi¬ 
cient to only consider quantum fluctuations of collective 
modes at low energies and long wave-lengths, e.g., den¬ 
sity and current fluctuations. Such an approach is the 
so-called “hydrodynamical approach” or effective quan¬ 
tum field theory for topological phases. The field theory 
realizations of SPT states belong to this approach. 

Historically, the “hydrodynamical approach” turns out 
to be extremely powerful to understand the underlying 
physics of topological phases. For example, the fractional 
quantum Hall effect (FQHE) can be understood by the 
Ginzburg-Landau Chern-Simons theory [32] or more sys¬ 
tematically by pure Chern-Simons theory [33-39]. Those 
bulk dynamic effective theories that capture the low en¬ 
ergies and long wave-length physics are also very useful 
to study phase transitions among different topological 


phases, e.g, phase transitions between FQHE at different 
filling fractions. Thus, the bulk dynamical Chern-Simons 
action approach to FQHE phases can be viewed as the 
Ginzburg-Landau action approach to symmetry breaking 
phases. Therefore, it is very natural to ask what is the 
“hydrodynamical approach” to SPT states. 

Very recently, Chern Simons/BF theories have been 
proposed [23, 40-45] as bulk dynamical effective ac¬ 
tions to describe 2D/3D bosonic SPT states protected 
by Abelian symmetry group (the so-called Abelian SPT 
states). Nevertheless, it has been pointed out[42] that 
the Abelian Chern Simons/BF theory approach is incom¬ 
plete. Thererfore, a much more general theoretical frame¬ 
work for bulk dynamical actions of SPT states is very de¬ 
sired. In this paper, we will focus on the mechanisms and 
bulk dynamical effective actions for bosonic SPT states 
with finite Abelian group symmetry within group coho¬ 
mology classification. We propose a class of new topolog¬ 
ical actions to characterize bosonic Abelian SPT states 
in arbitrary dimensions that are beyond Abelian Chern- 
Simons/BF theory. We will show that such a class of 
generalized topological actions serves as a complete de¬ 
scription for bosonic Abelian SPT states in ID and 2D. In 
3D, there are still some Abelian SPT states beyond the 
proposed bulk dynamical effective action; however, we 
believe that the basic principle and method developed in 
this paper are still applicable. We will leave these studies 
for future work. It is also worthwhile to mention that in 
a parallel work [46], a bulk dynamical effective action for 
Abelian SPT states beyond group cohomology classifica¬ 
tions is also proposed. In principle, the “hydrodynam¬ 
ical approach” can also be generalized into interacting 
fermionic systems. 


B. Summary of results 

1. A mechanism of SPT states 

Let us start by summarizing the mechanism that gen¬ 
erate SPT states at intuitive level. It is well-known that 
if we disorder a discrete-symmetry breaking state by con¬ 
densing domain walls, we can obtain a symmetry restored 
state. Our approach is basically analogous to this line 
of thinking, except that we generalize the approach by 
including additional multi-kink topological terms to the 
domain walls, see Fig.l. 

There are two ways to view the multi-kink topologi¬ 
cal terms: the space picture and the space-time picture. 
In the space picture, we create the symmetry-breaking 
domain walls and trap some charges (not fractionalized) 
of the remained unbroken symmetry at the intersecting 
points, then we proliferate and condense the domain walls 
to restore the broken symmetry. On the other hand, in 
the spacetime picture, we have an intersecting profile that 
contributes a nontrivial phase to the path integral (see 
Fig.l (b)), and we then disorder the symmetry breaking 
state with such nontrivial multi-kink topological term. 
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Dynamical gauge theory 

SPT invariants: 
Probed field theory 
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TABLE I: First column: the U k (l) non-linear er-model (NLuM) realization of the Zn 1 x Zn 2 x Zn 3 x ■ • ■ SPT states in 
the x < Xc disordered limit. The additional multi-kink topological term (bi-kink for 1+1D, tri-kink for 2+1D, quad-kink for 
3+1D, etc) are listed. The phase fluctuating term d^d 1 = d^dl + b £ contain a smooth piece dpOl and a singular piece b 1 ^. Here 
Cij... is a totally anti-symmetric tensor, with: C 12 = (^y , C 123 = (27r j 2 2! U 1234 = ( 2^331 * Plv . 

etc., with JV12... = gcd(A r i, A^ 2 ,...). Second column: the dynamical gauge theory realization of the Zn x x Zn 2 x Zjv 3 x • • • 
SPT states. The important global constraints on the fields are not specified, moreover we need to well-define the SPT path 
integral more than just the SPT Lagrangian; we will discuss this issues of path integral in depth in Sec. VII. Third column: the 
SPT invariants after integrating out the matter fields. Here the non-dynamical flat A 1 field describes the Z^-symmetry twist, 
which satishes <f A^dx^ = 0 mod 2n/Ni. The main result of our work is that the field theories in the first and the second 
columns are equivalent at low energies at the x < Xc disordered limit. We can derive their SPT invariants by integrating out 


the matter field. The SPT invariant is of the form: f d d x- 


-U1J2 


AjL^—^A^Ai 2 ■ ■ ■ Ai d given in [19]. 
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FIG. 1: (a) Disordering a [/(l)-symmetry breaking super¬ 

fluid with an action by condensing the vortices, e.g., tuning 
some coupling constant U to increase the charge repulsion. 
[47-49]. (b) Disordering a discrete-symmetry breaking state 
by condensing the domain walls. The gray region qualita¬ 
tively indicates the phase transition region, such as a critical 
point or a different phase, (c) In this work, we generalize the 
previous process by condensing domain walls with a multi¬ 
kink topological terms. We obtain nontrivial SPT states with 
SPT invariants listed in Table I. 


As we will show explicitly and quantitatively using field 
theories, both processes lead to a nontrivial SPT state. 

Using the above domain-wall condensation picture, we 
also obtain two kinds of field theory realization of the 
corresponding Z Nl x Z N2 x Z N3 x • • • SPT states (see 
Table I). The first one is a U k (l) non-linear a -model 
with a multi-kink topological term. The second one is a 
dynamical gauge theory that is beyond Abelian Chern- 
Simons/BF theory. Throughout the whole paper, we will 
implement the Euclidean spacetime approach with the 


Euclidean time tE = i t as the Minkowski time Wick- 
rotated by an imaginary i. We define the derivative do 
as dt E ■ We choose the Euclidean spacetime for the future 
convenience of the lattice regularization. 

In the first column of Table I, we list the U k { 1) non¬ 
linear cr-models with the nrulti-kink topological terms of 
the form 2 ^iJK...£ k ‘ v> ' 'd tl 6 I d u 9 J d\6 K ... with Cijk... a 
fully anti-symmetric tensor and d the spacetime dimen¬ 
sion. In the second column of Table I, we list the gauge 
theory realization of the same Z.v, x Zn 2 x Zn 3 x • • ■ 
SPT states. Our local field theories in the first and the 
second columns can produce the desired SPT invariants 
dictated by group cohomology [19] (after integrating out 
the dynamical fields). We list the SPT invariants in the 
third column of Table I. 


2. Field theory with anomalous U( 1) symmetry 

We stress that although the proposed U k ( 1) non-linear 
cr-model with the nrulti-kink topological terms formally 
has a U k (l) global symmetry 9j(x M ) —► 9i(x^) + Afj. 
Due to the presence of multi-kink topological terms, the 
U k ( 1) global symmetry is actually anomalous, i.e. can¬ 
not be realized by an on-site-synrmetry[ 12 ] in any lattice 
regularization of the field theories. Or more precisely, the 
C7 fc (l) non-linear cr-models have anomalous U k ( 1) sym¬ 
metry if the multi-kink topological terms are quantized 
as C 12 = 0 mod jX in 1+1D, C 123 = 0 mod ^ 2 -k) 2 2 \ 
2+1D, and C 1234 = 0 mod ^ 2 ^ 33 , in 3+1D. 

Here we use a 1 + ID example to explain the above 
statement (the higher dimensional cases can be under¬ 
stood in a similar way). Let us consider an ideal exper¬ 
iment by inserting a 27 t flux corresponding to the first 
1/(1) symmetry through a closed ID ring, the bi-kink 
topological term ^Cije^ v d ll 9 I d u 9 J will induce a charge 




































4 


27rC’i2 associate with the second U{ 1) symmetry. So if 
the 2ttCjj is not an integer, the U k ( 1) non-linear er-model 
does not even have the U k ( 1) symmetry at quantum level. 
When 2ttCij £ Z , the U 2 ( 1) symmetry is anomalous, 
since adding the flux of the first 17(1) can cause a non¬ 
conservation of the second 17(1). 

The above charge pumping phenomena via flux inser¬ 
tion can happen on a boundary of a 2 +ID system, where 
an integer charge is created in the bulk and the total 
U k (l) charges are conserved. 

However, the above charge pumping phenomena can¬ 
not happen in a strict 1 + ID system with on-site U k ( 1) 
symmetry. This is because the on-site 17 fc (1) symmetry 
is gaugable ( i.e. we can add 17(l)-flux without breaking 
the U k (l) symmetry). The presence of the charge pump¬ 
ing phenomena implies that, at quantum level, the U k (l) 
symmetry is broken by the 17(l)-flux, which in turn im¬ 
plies that the U k (l) symmetry is anomalous (or non-on- 
site). Or in other words, in a strict 1 + ID system with 
U( l) 2 on-site symmetry, C 12 must vanish. 

On the other hand, if the 2irCi2 = 0 mod N ^ 2 , 
the Z Nl x Zn 2 subgroup of the 17(1) 2 correspond to 
an anomaly-free symmetry (i.e. an on-site symmetry). 
This is because the 27r-flux of 17(1) induce a charge 
2 ttC \2 = yvW x integer, which is essentially trivial since 
Zn 2 charge is only conserved mod N 2 . Therefore, the 
Zm 1 x Zpf 2 subgroup of the 17(1) 2 is no anomalous. 
The 17 2 (1) non-linear cr-model describe a system with 
Zn x x Zn 2 on-site symmetry, if 2 nCi 2 = 0 mod yyy 2 -- 


Similarly, the U k (l) non-linear cr-model have a Z jv, x 
Zn 2 x Zj\r 3 x • • • on-site symmetry only when proper 
quantized values is assigned for Cijk...- For example, 
in 1 + ID, 2 + ID and 3 + ID. we require that C 12 = 
1 .Mm Ci23 = and Cl23 4 = 


(2tt) 

1 

(27r) 3 3! 


JV 12 ’ ^ ~ (27r) 2 2! N 123 

NlN2 N^:' i P1V > where PhPihPm e z. 


C. Organization of the paper 

The rest of the paper is organized as follows: In Section 
II, we briefly review how to use SPT invariants to define 
SPT states. In Section III, we propose a bulk dynamical 
effective action to describe 1+1D bosonic Abelian SPT 
states and use it to derive the corresponding SPT invari¬ 
ants. In Section III, we briefly review the Chern-Simons 
action approach for 2+1D bosonic Abelian SPT states 
and discuss its limitation. In Section IV, we compute 
the SPT invariants for 2+1D Zn x x Zn 2 x Zjm 3 SPT state 
and propose a bulk dynamical effective action to describe 
such 2+1D SPT states. In Section VI, we generalize our 
results to 3+1D bosonic Abelian SPT states and pro¬ 
pose a bulk dynamical action beyond BF theory. In Sec¬ 
tion VII, edge theories for Abelian SPT states beyond 
Clrern-Simons/BF actions are discussed via a standard 
dimension reduction scheme. Finally, there are conclu¬ 
sion remarks and discussions for future directions. 

In Appendix A, we review the derivation of disordering 


the superfluid state to the Mott insulator, see the pioneer 
work[47-49] and Ref. [50, 51]. In Appendix B, we provide 
an explicit calculation of an effective bulk action of SPT 
state. In Appendix VII, we verify that the partition func¬ 
tion with the proposed SPT action has the GSD=1. In 
Appendix C, we provide some words of caution by com¬ 
paring our effective action of SPT state to topological 
gauge theories with non-semi-simple Lie algebra. In Ap¬ 
pendix D, we compute the edge mode GSD by counting 
the degenerate zero modes. 


II. A REVIEW OF SPT STATES DEFINED BY 
SPT INVARIANTS 

It has been shown that SPT states (within group co¬ 
homology or beyond group cohomology classifications) 
can be probed or even defined through the so-called SPT 
invariants [21, 52] that may completely characterize dif¬ 
ferent SPT states. In this section, we will review and 
discuss such a point of view. 


A. Universal wavefunction overlap: a complete 

SPT invariant for SPT orders 

We start from reviewing the results of the SPT invari¬ 
ants in Ref. 52, using 2+1D systems as examples. It 
was conjectured that the degenerate ground states |4' a ), 
a = 1, 2, • • •, of a 2+1D topological phase on a torus have 
the following properties: [53] 

S af B e~f sL2+ °^ = <* a |£|*0> 

T a p e-A^+oCL- 1 ) = (* Q |f |^) (i) 

where S is the 90° rotation operation (x,y) —» (—y,x) 
and T is the Dehn twist rotation operation (x, y) -A 
[x+y, y). It was conjectured that while the complex num¬ 
bers fs and /t are not universal, the complex matricies 
S a p and T a p are universal. S a p and T a p can change only 
via phase transitions. Thus we can use them to charac¬ 
terize different topological orders. In fact, we believe that 
S a /3 and T a/3 completely define 2+1D topological ordered 
phases with gappable edges. 

Can we use the similar idea to completely define 2+1D 
SPT order? The wavefunction overlap for SPT state also 
has the following universal structure 

S e - fsL2+ol ' L ~ 1) = ('FoISj'Fo) 

T e ~ fTL2+o(L ~ 1) = ('Folij'Fo) (2) 

where the lxl unitary matrices S and T are universal. 
In fact S = T = 1, due to the trivial bulk excitations in 
SPT state. Thus S and T are trivial and could not be 
used to distinguish different SPT states. 

To obtain a non-trivial wavefunction overlap, we intro¬ 
duce symmetry twist: a symmetric transformation gen¬ 
erated by h £ G within the region R. The Hamiltonian 
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FIG. 2: (a) Symmetry twist along the boundary dR is gen¬ 

erated by the symmetry transformation that act only within 
R. (b) The symmetry twist h x ,h y on torus gives rise to the 
twisted ground state 



FIG. 3: S'-move is 90° rotation. 



FIG. 5: A closed orbit in the (h x ,h y ) space. 

O 

h x 

FIG. 6: Two closed orbits in h x space. 

The S-, T-, and 17-moves shift ( h x ,h y ) (h' x ,h' y ) 


is not invariant under such a local symmetry transforma¬ 
tion (see Fig. 2(a)): 

H = J2 H ijk ->■ H h = J2 H W + E H ijk (3) 

in R, R on dR 

where H l]k acts on sites i,j, k and H^ k is on the bound¬ 
ary of R , dR, if the sites i,j, k are not all on one side of 
dR. We call ^) on gR H^ k the /i-symmetry twist. 

Note that H and Hh have the same energy spectrum. 
So the symmetry twist costs no energy. Let | ^(h x ,h )) 
be the ground state of Hh It h v on a torus with symmetry 
twists h x ,hy in x- and y-directions. |’F(fc, x ,/i !; )) simulates 
the degenerate ground states for topologically ordered 
phases. We can use )) to construct S,T matrices 

that characterize the SPT order (see Fig. 3 and Fig. 4): 
S move: {^ {h -\ h JS\^ {hx>hy) )=S hxthy e ~fsL 2 +o(L - 1 ) 

T move: (* { h x ,h y hjT\* {hx ,h y )) = T K ,h y e~^ La +°^ 
U move: {^ {hthxh -\ hthyK ^\U(h t )\^(h x ,h y )) = U h ^ hy (h t ) 

Note that in addition to the S- and T-moves, the SPT 
invariants also contain U- move generated by the global 
symmetry transformation h t £ G. 


S : ( h x , hy) > (h x , hy) — (h y , h x ); 

T . {h x ,h y ) ^ (.h x ,h y ) {h x , h y h x ), 

U(ht ) : {hx, hy) ->■ {h'x, h'y) = (hthxhi^hthyhi 1 ). (4) 

When ( h' x ,h ') / (h x ,h y ), the complex phases 

Sh x ,h y ,T hx ,h v ,U h:D 'hy{ht) are not well defined, since 
they depend on the choices of the phases of 
\^{h x ,h y )) and | ^{h> x ,h' y ))- However, the product 
of S hiC ' hy ,T h:C ' hy ,U hxihy (h t ) around a closed orbit 
(h x ,h y ) —> (h x ,hy) —>■•••—)• ( h x ,h y ) is universal (see 
Fig. 5). We believe that those products for various closed 
orbits completely characterize the 2+1D SPT states. 

For example, N T-moves always form a closed orbit 
for Abelian Zjy = {h = 0, • • • ,77 — 1} group. For 2+1D 
Z N SPT state labeled by k £ H 3 [Z N , t/(l)] = Z N , it has 
one SPT invariant: 


T, 


h x 


■■■Tv. 1,2 u Tv. h h Th h = e 27ri (^- 1 ) 2fc / iv 

h h x h z ,h y - 1 - h x h y ,h y h x ,h y — c 5 

( 5 ) 


h X 5 hy £ 


N • 


Such an SPT invariant completely characterizes the 
2+1D Z N SPT state. 



(a) (b) (c) 


FIG. 4: T-move is the Dehn twist followed by a symmetry 
transformation h x in the shaded area. 


B. Universal wavefunction overlap in 1 + 1D 

In 1+1D, the SPT invariants are very simple. We only 
have the U-move: {hthxh -^\U[h t )\^ {hx) ) = U hx {h t ), 

which generates the shift h x —> hthxh ^ 1 . Similar to the 
2+1D cases, the product of Uh x (h t ) around a closed orbit 
is well defined and universal (see Fig 6). In particular, 
for Abelian symmetry group, Uh x {ht) itself is universal. 
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FIG. 7: Space-time S 1 x S 1 with two symmetry twists in x,t 
directions. 


ttT^TTT 

Jz ^z 

FIG. 8: Two kinds of domain walls with the same energy, but 
different Zf-charges, 0 (mod 2) and 1 (mod 2) respectively on 
a lattice. Eq.(14)’s H^° p is a hopping operator for the first 
kind of domain wall. Eq.(15)’s is a hopping operator for 
the second kind of domain wall. 



III. A 1+1D Z Nl x Zjv, SPT STATE AND ITS 
BI-KINK BULK DYNAMICAL ACTION 


A. A simple example 


Now let us apply the results obtained in the last Sec¬ 
tion to a 1+1D Zj\q x Zm 2 bosonic SPT state, which is 
classified by 

H 2 [Z Nl x Z N2 ,U{ 1)] = Z Nl2 = {0,1, • • • , N 12 - 1} (6) 


where iVi 2 = gcd(7Vi, iV 2 ). We consider an SPT state la¬ 
beled by k £ Z at 12 . 

The group elements of Zjvj x Zat 2 are labeled by 
h = (h 1 ^ 2 ), h 1 € Z Nl , h 2 £ Zjv 2 . The universal wave- 
function overlap (the SPT invariant Uh x (ht)) is 




= U h i h i{h\,hf) = e 


ik wk( h l h t- h l h l) 


( 7 ) 


which can also be viewed as the fixed-point partition 
function on space-time T 2 = S 1 x S' 1 with symmetry 
twists in x,t directions (see Fig 7): 


^fixed-point = U h l tK (hl,h 2 ) = (8) 


The twisted ground state carries a non-trivial quantum 
number e ,lr ^ Si with e I7r ^ Si = — 1. Such a dependence 
of the ground state quantum number e™ ^ Si on the 
twist is the 1+1D SPT invariant discussed above. 

The above SPT invariant also suggests a mechanism for 
the 1+1D Zn x x Zn 2 SPT state. We notice that the SPT 
invariant implies a symmetry twist of Z^q that carries a 
“charge” of Z N2 . Since the symmetry twist of Z Nl is the 
domain wall of Z^q in a% symmetry breaking state, we 
may (1) start with a Z^ l symmetry breaking state, (2) 
bind k Zjv 2 -charge to the domain wall of Z^q, and (3) 
restore the Z^q symmetry by proliferating the domain 
walls. In this way, we obtain a 1+1D Z^q x Zn 2 SPT 
state labeled by k £ / H 2 [Zn 1 x Zn 2 , t/(l)]. 

For example, let us consider a ID Zf x Zf spin-1 chain 
with symmetry 

Z% : (| t*>, |0x), | -U» -t (| tx), -|0 X ), | lx)) 

Zl : (I tz>, |0,), | -A (| fz), -|0,), | 4-z))- (12) 

The following Hamiltonian has the Zf x Zf symmetry 

H°z 2X z 2 =T,- J * S t S U i ( 13 ) 

i 


Both wavefunction overlap and partition function pic¬ 
tures imply the following physical meaning for the above 
SPT invariant: a symmetry twist of Zj\q carries Zjy 2 - 
charge k. 

^{hl,hi)={ifi)\U[h] = 0, h 2 = l)| , Fp l i,/ l 2) = ( 1 , 0 )) = e lfejv i2, 

( 9 ) 

Let us discuss a concrete example for the above 1+1D 
SPT invariant. We consider a spin-1 chain with the spin- 
rotation symmetry Z 2 x Z 2 = D 2 = 180° in S x . S z . The 
Hamiltonian on a ring is given by (untwisted case) 

L — l 

H D2 = ^(J X S*S? +1 + j y s*sv +1 + J z s z s z +1 ) 

2=1 

+ J x S x L Sf + J y S v L S\ + J z S z L Sl (10) 

where J x = J y = J z > 0. The ground state carries a 
trivial quantum number e l ' K ^ Si with e l7T ^ Si = 1. 

If we add a twist by e I7r ^ s *, the Hamiltonian becomes 

L—X 

rrtwist _ \ ' / T QX QX I T Qy QV i T Q z Q z \ 

n D 2 ~ 0 i+l J V D i J i+1 ' J z J i J i+1) 

j=l 

+ J x S x L Sl - J y S v L S\ - J z S z L Si (11) 


but its ground state breaks the Zf symmetry. Such a 
symmetry breaking state has two kinds of domain walls 
which happen to have the same energy, but different Zf- 
charges. The two kinds of domain walls, shown in Fig. 8, 
have different hopping operators: 

h\ op =£i-f(OT + m 2 ) 

= Ei-K((Sf) 2 -(sy) 2 ), (14) 

tf 2 hop =-^^(^+1 + ^5^) 

= £ iJxvi-Sf^ + SVS?^). (15) 

Here we used the fact that S+ = Sf + i Sf and S~ = Sf — 
iSf. It is straightforward to see the (S+) 2 operator hops 
the first kind of domain wall of Fig. 8 in one direction, 
while the ( S ~) 2 operator hops the first kind of domain 
wall of Fig. 8 in the opposite direction. On the other hand, 
the S +operator hops the second kind of domain wall 
of Fig. 8 in one direction, while the S~ S~ +1 operator hops 
the second kind of domain wall of Fig. 8 in the opposite 
direction. 

Adding a strong enough hopping operator can make 
a domain wall subject to a negative energy cost, which 


7 


restores the symmetry by proliferating the domain 
walls. We find that Hz 2 xz 2 + Hi° p leads to a trivial 
SPT state, while Hz 2 xz 2 + ^° P leads to a non-trivial 
Z -2 x Z 2 SPT state. Via a unitary transformation, the 
Hamiltonian H ZoxZo + H 2 ° p is equivalent to the Hamil¬ 
tonian of Eq.(10) discussed above, as the Haldane phase 
of a spin-1 anti-ferromagnetic Heisenberg chain. 


B. Bi-kink topological term NLcrM and dynamic 
gauge theory 

The underlying physics of the above 1+1D Z^ 1 x Z^ 2 
SPT state can also be captured by the following Higgs 
action with a bi-kink topological term: 

£ bi -kink = | {d ll e I ) 2 + 1 -c IJ e^d li e I d v e J 

m+w (i 6 ) 

+ -C' / 'A A ‘ ! '(<9 M 0j + + K) + ^Maxwell’ 

where 1 = 1,2 and the structure constant Cij is totally 
antisymmetric with Cu = — Cjj. We assume Einstein 
summations for repeated indices throughout the whole 
paper. The quantum phase fluctuation can be captured 
by a real scalar compact field 9 1 = 9f + 9f with a smooth 
piece and a singular piece 9\ and 9\. To achieve the 
disordered insulator state, we can condense the vortex, 
namely strongly disorder the superfluid coherent phase. 
We will write d^9l +d ll: 9{, = d^9l +6^. The captures 
the smooth piece d and the additional b 1 captures 
the singular piece d^9f. We note that the real scalar 
fields 9 J S can be viewed as the phase fluctuations of x7y f 
symmetry in a, Zn 1 x Zn 2 symmetry breaking phase while 
vector fields b 1 (with T^ axwell the corresponding Maxwell 
term) describe the proliferations of domain walls, which 
restore the Zj^ 1 x Zy 2 symmetry. Such a Higgs action 
with a bi-kink topological term will enforce a Z^ 1 domain 
wall that carries a “charge” of Zn 2 , and vice versa. It is 
clear that the bi-kink topological term is just a boundary 
term in the absence of gauge fields b^. In the following 
we will show that such a bulk action with the bi-kink 
topological term indeed describes the Zm 1 x Zn 2 SPT 
physics in 1+1D. 

After dropping the total derivative term, we can 
rewrite the above action as: 

£bi-kink = f(9X+^) 2 ( 17 ) 

+ 2 C ,/j£/iy ( — + b^b^) + ^Maxwell! 

Next, we introduce the Hubbard-Stratonovich fields jj 
to decouple the quadratic term as 

£bi-kink = + 'tfji 

+ -^C IJ ^ V { — ^ I s 3^b J v + b 1 ^) + /^Maxwell) 


Integrating out the smooth fluctuations 9f leads to the 
following constraint: 

d^ + C IJ e^b J v ) = 0. (18) 

The above constraints can be solved by: 

f = ^e^dvd 1 - C IJ £^b J „. (19) 

Z7T 

where a 1 do not need to be globally defined. To disorder 
the U( 1) phase, we take \ <C Xci we can drop out the 
5^0/) 2 term as well as the Maxwell term of gauge helds 
b 1 thanks to their RG irrelevancy [51]. We end up with 
an effective topological action: 

Aop = ±e' t %d v a I + ^C IJ e»%b J „, (20) 

The gauge transformation of b 1 in the above action will 
induce a shift on the scalar fields a 1 : 

a 1 —■► a 1 + 2nC IJ g J -, b 1 ^ —» b 1 ^ + d^g 1 . (21) 

The above functions do not necessarily need to be globally 
defined. In fact, the compactness condition of a 1 and b 1 
implies the closed loop or the closed surface integral has 
the constraints: 


/ 


da/( 27 r) € Z, 


dbj/(2n) G Z. 


( 22 ) 


In Sec. VII, we will derive the same constraints in the 
path integral level, from the constraints of 17(1) charge 
and the vortex number on a closed-surface. 

Now, let us compute the quantization condition for 
the coefficients Cjj. We note that the average of 9 1 = 
9f + 9 x is quantized as 27r/JV/X integer. In the disordered 
phase which restores the Z^ 1 x Z^ 2 symmetry, 9 1 's have 
many fluctuating kinks. Let us consider a configuration 
where 9 1 has a kink A 9 1 = 2nki/Ni on the t axis and 
9 2 has a kink A 9 2 = 2irk 2 /N 2 on the x axis. For such 
a configuration, the action from the bi-kink topological 
term is given by 


S = 


dxdt -Cue^d^b 1 d v 9 J 


= 87r 2 i C 


12 


k\k 2 

N X N2 


(23) 


This means that the 9 1 kink carries a iyv 2 -charge 
27rCi2^j- mod N 2 . Since k\ = 0 ~ ki = Ni, C \2 must be 
quantized: 


2nCi2 = 0 mod N 2 , 2ttCi 2 = 0 mod Ni. (24) 


Thus 


_ Pn N\N 2 
2n N 12 


Pn — 0, • • • , N 12 — 1 


(25) 




where N i2 = gcd(iVi,iV 2 ). Also we note that C 12 has 
only N 12 distinct quantized values, corresponding to Ni 2 
distinct charge assignments. 

The above argument for the quantization condition of 
C\ 2 due to global Zn 1 x Zn 2 symmetry can also be de¬ 
rived in a rigorous way by adding a coupling term to 
external background gauge field A 1 : 

^coupling = ^s^A^a 1 , (26) 

As the physical meanings of A 1 and A 2 are Z^ x and 
Zjy 2 symmetry twists, A 1 must be a flat connection with 
d A 1 = 0 and j> A 1 = 2nnj/Nj. On the other hand, since 
f dx dt£ COU piing must be invariant under gauge transfor¬ 
mation Eq.(21), C 12 can not take arbitrary value. A 
shot calculation suggests the same quantization condi¬ 
tion Eq.(24). 

In Sec. VII, we will define a rigorous SPT internal gauge 
theory path integral, and we confirm that the GSD of our 
theory is unique on a closed manifold, GSD=l,in agree¬ 
ment with SPT state. We will also derive the SPT in¬ 
variant in Ref. [19] by coupling the internal gauge theory 
to semi-classical probed field A. In this way, it becomes 
manifested that C 12 can only take N- t 2 distinguishable 
value derived in Eq.(24). In the following, we generalize 
the above results to higher dimensions. 

IV. A REVIEW OF CHERN-SIMONS ACTION 
APPROACH TO 2 + ID ABELIAN SPT STATES 

In this section, we will start with a brief review on 
the Chern-Simons action approach for 2+1D Abelian 
SPT states. Then we explain the physical meaning of 
the Chern-Simons action approach and discuss its limi¬ 
tations. 

It is well known that a vortex condensation can turn 
a boson superfluid into a trivial bosonic insulator. A 
bosonic U( 1) SPT state is also a bosonic insulator, but 
a non-trivial one. It turns out that a condensation of 
vortex-charge bound state can turn a boson superfluid 
into a non-trivial U( 1) SPT state. 

To show this, let us consider a boson superfluid for 
one specie of bosons, which can be described by an XY 
model: 

C X Y=^(d^0) 2 , (27) 

If the vortex of the boson condenses, 6 in the XY model 
is no longer a smooth function of space-time. We can 
introduce the singular part by replacing d^O by d^Og + b^, 
where the field strength of gauge field corresponds to 
the vortex current density = ^£^ vX d v bx. 

The charge of gauge field is the number of vor¬ 
tices minus the number of anti-vortices and is quantized. 
In the vortex condensed phase, the phase fluctuation of 
the vortex condensate can be described by another XY 
model, which is dual to the Maxwell term of the gauge 


field 6 m . Now the boson superfluid is described by the 
following Lagrangian 

Auggs = \[(d^e s + A ) 2 + 4^2 vn (28) 

where — d v b M and we have normalized with 

v = 1,X = !■ 

We can introduce a Hubbard-Stratonovich field j fI to 
decouple the quadratic term as 

Auggs = ~ + i b,f + ^dis¬ 

integrating out the 0 S field results in a constraint = 
0. From this constraint, we can write = -^-£ ,lvX d v a\. 
The charge of a^ is equal to the boson number and is 
quantized. With these results, the path integral becomes 

£bf = i^ A M^ + i[d^ + ^H,(29) 

where F flv = — d^a^. Note that the boson cur¬ 

rent j b = ^£ >ivX d v bx, while the vortex current j” = 
A; £^ x d u a\. 

The above can be generalized to the case with fc-species 
of bosons with C7 fc (l) symmetry. The bosonic insulator 
induced by the vortex condensation is described by the 
following Chern-Simons action: 

C C s = ^ vX K 0 IJ a I lx d v a{-, I = 1,2, • ■ • , 2k (30) 
with 

K ° l1 = ( 1 0 ) ® Ifexfe (31) 

where a 2 ^ ~ a^ and a ^ fc_1 ~ b /i . Since | det[A']| = 1, the 
above Chern-Simons action has a unique ground state 1 
on any closed manifold. The chiral central charge for the 
edge states is given by the signature of K which is zero. 
So the bosonic insulator has a trivial topological order. 

However, the bosonic insulator may have a non-trivial 
U k ( 1) SPT order. To see this, we turn on the external 
U k ( 1) gauge field A“ to reveal the U k ( 1) symmetry of 
the theory: 

Coupling = a = 1, 2, ■ ■ • , k (32) 

Here q a are integer-value charge vectors. is the 

A“-charge carried by the T h -species of bosons, and qjf is 
the A“-charge carried by the vortex of the Z th -species of 
bosons. We see that charge vectors q a reveal the infor¬ 
mation on what kinds of vortex-charge bound states are 
condensing to produce the bosonic insulator. Different 
vortex-charge bound states (i.e. different charge vectors) 
will lead to different U k ( 1) SPT orders. 

The full theory is given by £ = £ C S + ^coupling, Af¬ 
ter integrating out internal gauge fields (the matter 
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fields), we obtain an effective theory for the external fields 
A a : 

As = --^s^ x AXK°uqjd v Al (33) 

By considering equivalent class of response K matrix 
K a p = qaK°uqp, we can “classify” 2+1D U k ( 1) SPT 
states described by the Chern-Simons theory Eq.(30). 
We can also break the U k (l) symmetry down to Zn x x 
• • • x Z]y k symmetry and obtain a “classification” of 
Zn 1 x • • • x SPT states in 2+1D. Since Zu a group 
can always be embedded into U k ( 1) group, it is not a sur¬ 
prise that the Z^, x ■ ■ ■ x Z^i k SPT state can be described 
by the same Chern-Simons action. 

However, since H 3 [Z Nl x • • • x Z Nk ,U( 1)] = 

I-Nij ®i<j<k Z N ijk { N ijk = gcd (Ni,Nj,N k )), the above 
classification turns out to be incomplete and it can only 
describe a subclass of Abelian SPT states labeled by 
®iZ-Ni ®i<j Zjv y , namely, the type I and type II SPT 
phases. In the following, we will develop an effective 
field theory description for type-III SPT order in 2+1D, 
which is labeled by ® i<j<k Z N .. k . 

V. A 2+1D Zn i x Zn 2 x Z N3 SPT STATE AND ITS 
TRI-KINK BULK DYNAMICAL ACTION 

A. A 2+1D Zn x x Zn 2 x Zn 3 SPT state 

Without loss of generality, it is sufficient to discuss 
a 2+1D Zn 1 x Z N2 x Z Na bosonic SPT state, which is 
classified by 

H 3 [Z Nl xZ N2 x Z Ns ,U(l)) (34) 

= Zjvq x Zjv 2 © Zjv 3 © Z Nl2 © Zjv 23 © Zat 13 © Z n 123 

We consider a type III SPT state labeled by ft £ Zjv 123 

The group elemenets of Zn 1 x Zn 2 x Zn 3 are labeled 
by ft = (ft 1 , ft 2 , ft 3 ), ft 1 £ Zjvq, ft 2 £ Zjv 2 , h 3 £ Zjv 3 . The 
SPT invariant U kx h (ft t ) for the above SPT state is the 
fixed-point partition function on space-time T 3 = (S' 1 ) 3 
with symmetry twists in x, y , t directions: 

^fixed-point = U hx , hy (h t ) = (35) 

The physical meaning of the SPT invariant is the follow¬ 
ing: Consider the ground state of the Hamiltonian with 
symmetry twists in Zpj 1 and Zn 2 , the intersection of the 
symmetry twist in Z^ 1 and the symmetry twist in Zk 2 
carries Zk 3 -charge k. 

The above SPT invariant also allows us to calculate the 
dimension reduction of the 2+1D SPT state to a 1+1D 
SPT state: We view the space-time as T 3 = T 2 t x S' 1 , 
and put Zn 3 symmetry twist (ft. 1 , ft)), ft 3 ) = (0,0,1) in the 
small circle S 1 . The 2+ID partition function reduces to 
a 1+1D partition function 

^fixed-point = e ik2 ^ h * h *- h * h ^ (36) 


which is the SPT invariant of a 1+1D SPT state. We find 
that the resulting 1+1D SPT state is the one labeled by 
k £ H 2 [Z Nl x Z N2 ,U(1)] = Z Nl2 . The boundary of such 
a 1+1D SPT state carries degenerated states that form a 
projective representation of Z.v, x Zn 2 . This leads to an 
experimental probe of the Zn x x Zjy 2 x Zn 3 SPT state: 
a Zpf 3 “vortex” (end of Z N:i symmetry twist) carries de¬ 
generated states that form a projective representation of 
Z N i x Zm 2 - 

The result of the above dimension reduction can also 
be viewed as each z n 3 twist (which is a ID curve in 
2D space) carries a 1+1D Z^ x Zm 2 SPT state labeled 
by k £ H 2 [Zm 1 x Zjv 2 ,U(1)]. This picture leads to an¬ 
other mechanism for the 2+1D Zk x x Z^ 2 x Zjv 3 SPT 
state: (1) start with a Z^ 3 symmetry breaking state, (2) 
bind a 1+1D Zk 1 x Zk 2 SPT state to the domain wall 
of Zn 3 , and (3) restore the Zu 3 symmetry by prolifer¬ 
ating the domain walls. In this way, we obtain a 2+1D 
Zn 1 x Zn 2 x Zn 3 SPT state labeled by ft £ Zjv 123 . 

The 2+ID SPT invariant Eq. (35) on space-time 
T 3 = (S' 1 ) 3 can also be expressed as a topological term 
of probe fields Aj: 


/vtwist ( / T'3 
^fixed-point 


= e 


ipin 


N 1 N 2 N 3 
(2 7r£iV 123 


f AiAA 2 AA 3 


d Aj = 0, 
(37) 


with an integer pui- Again, since Aj describes symme¬ 
try twists on the boundary, it must be flat connection 
with dA/ = 0. f Ai A A 2 A A 3 is also gauge invariant 
if d Aj = 0. The field theory representation of the SPT 
invariants Eq. (37), should be valid for any space-time 
topologies. In the following we will show how to derive 
such a topological response from a bulk dynamical effec¬ 
tive action. 


B. Tri-kink topological term NLcrM 

To describe the so-called type-III Z Nl xZ N2 xZ Na SPT 
orders in 2 + ID, we consider the following effective ac¬ 
tion for three species of bosons with vortex condensation. 
The action contains a new tri-kink topological term - the 
Cijk -term (the following is a generalization of Eq. (28)): 

Ari-kink = l^e 1 ) 2 + ±c I j K e*'%o I d v o J d x e K 

— 2 (A^s + ^m) 2 + -^Maxwell (38) 

+ \c IJK e^\d,ei + bl)(d^ + K)(d x ef + ftf) 

where I = 1, 2, 3 and the structure constant Cijk is to¬ 
tally antisymmetric with Cijk = —Cjik = —Cikj- It 
is clear that the tri-kink topological term is just a bound¬ 
ary term in the absence of gauge fields b 1 . 

To understand the physical meaning of the tri-kink 
topological term, we first note that the type-III SPT or¬ 
ders in 2+ID only exist for a finite group Zpj 1 x Zjy 2 xZk 3 - 
So we need to break the t/(l) 3 symmetry down to 
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Zn-l x Zn 2 x Zn 3 symmetry. The average of 9 1 is quan¬ 
tized as 2ir/NiX integer. In the disordered phase which 
restores the x Z N2 x Zj^ 3 symmetry, 9 T, s have many 
fluctuating kinks along space-time surfaces. Let us con¬ 
sider a configuration in the space-time where 9 1 has a 
kink A 9 1 = 2-Kk\/N\ on the y-t plane, 9 2 has a kink 
A 8 2 = 2Trk 2 /N 2 on the t-x plane, and 9 3 has a kink 
A 9 3 = 2nk 3 /N;i on the x-y plane. For such a configu¬ 
ration (b^ = 0 ), the action from the tri-kink topological 
term is given by 


S= dxdydt l -C I j K £^ x d ll 8 I d„9 J dx9 K 


= 167r d i C 


123 


k\k 2 k 3 
N\N 2 N 3 ■ 


(39) 


This means that the intersection of the kinks in 9 1 and 
9 2 carries a Zn 3 - charge 87r 2 C l i 2 3 ' mod N 3 . Since 

fci = 0 ~ k\ = Ni, C 123 must be quantized: 


u fci 

87 r 2 Ci 23 —— = 0 mod N 3 , 87 t 2 C'i 2 3 — = 0 mod N 3 . 
N 2 N i 


(40) 


Thus 


Pin N 3 N 2 N 3 
(2tt) 2 2! A 123 


Pill = 0, • • • , N 123 — 1 

(41) 


where N 323 = gcd(N 3 ,N 2 ,N 3 ). Also we note that C 123 
has only IV 123 distinct quantized values, corresponding to 
A 123 distinct charge assignments. 

Now the physical meaning of the tri-kink topological 
term is clear: It is well known that the fluctuations of the 
kinks will turn a Zj^ 1 xZ^ 2 x Zjy 3 symmetry-breaking 
state into a Zjyj x Z^ 2 x Zjy 3 symmetric state with a 
trivial SPT order. However, if we bound a Z^-charge 
to the intersection of the kinks in 9 1 and 6 2 etc , the 
resulting Zn x x Zjy 2 x Z^r 3 symmetric state will have a 
non-trivial SPT order, as we will show below. In this way, 
we can produce N 123 distinct type-III Zn x x Z-^ 2 x Z^ 3 
SPT orders, consistent with the group cohomology result. 

By integrating out the smooth fluctuations 8 1 and in¬ 
troducing auxiliary gauge fields a x and A^, we can derive 
the following bulk dynamical action: 


Ari- kink = + \c IJK e» vX { A^Af + (&£ - A£)(tf - A J v ){b* - Af)) + - A ^) 2 + Tl^weii- (42) 


The derivation from Eq.(38) to Eq.(42) is preserved in Appendix B with details. Interestingly, the field strength 
of gauge field a^ is formally akin to a non-Abelian gauge field and its infinitesimal gauge transformation should be 
modified as: 


+ - IttCuk ( 9 J A* + ^g J d ll g 


K 


bji —■> b li + 9/j.g 1 ; A^ — > A^ + d^g 1 . 


(43) 


C. Saddle point approximation and internal gauge 
theory 


If we assume the held b 1 has a weak fluctuation, we 
can apply the saddle point approximation for b^. The 
saddle point equation reads: 


and with the gauge redundancy given by: 

K -> bi + d^g 1 ; 

-> + d^f 1 - A-kCuk (g J b* + ^g J d^g K ^j . (46) 

We also have the global constraints: 


C IJK e^ x {bi - A i){b* - Af) + (&£ - A£)+ 

higher order terms = 0, (44) 

clearly b 1 = A^ is a stable saddle point. Since the A 
held is a Lagrangian multiplier and b is a more-restricted 
17(1) held, we should replace A by b. At the level of this 
approximation, we can simplify the bulk effective action 
by: 


da/(2n) G Z, d&//(27r) £ Z (47) 

Similar to the 1 + ID case, there is a rigorous way to 
compute the quantization of coefficients C 123 protected 
by global ■Z’/v, x Zn 2 x Zjy 3 symmetry. Let us add a 
coupling term to the external gauge held A 1 . 

71 coupling = iAtf = 7 ^£^ x A I ll d iy a[, (48) 



£eff = 


1£ 


His\ 


2n 


b ^a x 


+ ^e^bibib f 


A 7 


(45) 


Again, A 1 are Zn i symmetry twists, thus A 1 must be 
a hat connection with d A 1 = 0 and f A 1 = 2 tttij /Nj . 


3 
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Similar to the 1 + 1D case, since f dxdydt£ coup u ng must 
be invariant under gauge transformation Eq.(46), C 123 
can not take arbitrary value, and a short calculation gives 
rise to exactly the same condition Eq.(40). 

It turns out that the above gauge transformation cor¬ 
responds to a non-semisimple Lie algebra of symmetry. 
We will discuss a generic class of such Lie algebra, called 
the symmetric-self dual Lie algebra in Appendix C. In 
Sec. VII, we will define a rigorous SPT internal gauge the¬ 
ory path integral, and we confirm that the GSD of our 
theory is unique on a closed manifold, GSD=1, just like 
the SPT state. We will also derive the SPT invariant 
by coupling the internal gauge theory to semi-classical 
probed field A claimed in Ref. [19], which suggests that 
Eq.(41) indeed gives rise to V 12 3 distinguishable SPT 
phases. 

VI. A 3+ID GENERALIZATION 

The above tri-kink topological term can be general¬ 
ized into higher dimensions as well, such as a quad-kink 


topological action in 3+1D: 


£ q -ki„k = * * (9X + K) 2 + \c IJKL e^{dX + bl) 

■{d„ei + btWtf? + b«)(d a eL + b L a ). (49) 


The quantization condition on Cijkl can be worked 
out in a similar way, and finally we obtain C 1234 = 
where piv is an integer on p w = 

0 , ■ • • , V1234 — 1 - 

For example, in 3 + ID, we can use the following 
quartic-kink term to describe the so-called type-IV SPT 
state. Parallel to our previous derivation in Sec.VB, we 
can derive the SPT bulk dynamical action: 


£ q -kink = 

and its gauge transformation: 




^){b L p -\ L p )) + ~(b 


\I \2 1 rb 
/x [iJ ' ^ Maxwell 


(50) 


a i.v Jl 


dvff l +‘2dTrC I jKLg J ^^u + ---', b 1 ^ —> b 1 ^ + d^g 1 + ...; -> + d p g J + .... (51) 


Here ^Maxwell t erms contain non-topological Maxwell 
term. If we further apply the saddle-point approxima¬ 
tion, we obtain: 

£eff = b L p . (52) 

The corresponding infinitesimal gauge transformation(we 
only keep the leading order term here and use ... to 
represent higher order terms) of arbitrary functions / 
and g reads: 

^ a U + d rfl - d v fl+2ATTC I j KL g J b 1 ^b^+ ■ ■., 
b^ t b p T d fl g I + .... (53) 

Here g 1 and b 1 are globally defined, but f 1 is not globally 
defined. The analogous global constraint can be derived: 

da/(27r) e Z, $ dh//(27r) £ Z. (54) 

VII. PARTITION FUNCTION, GSD AND SPT 
INVARIANTS COMPUTED FROM THE SPT 
INTERNAL GAUGE THEORY 


Here we will analytically show the path integral defi¬ 
nition of internal gauge field theory, Eq.(20) for 1+1D, 
Eq.(45) for 2+1D, Eq.(52) for 3+1D. In particular, we 
will show three key issues: 

• (f) Define the partition function Z using field theory 
path integral. 

• (ii) Derive the SPT invariants of semi-classical flat 
probed field theory in Ref. [19] by coupling the SPT in¬ 
ternal gauge theory to probed fields A. 

• (Hi) The internal gauge theory on any compact closed 
spatial manifold has a unique ground state, namely 
GSD=1. This means that the absolute value of the phase 
space volume ratio between the case with topological 
term and the case without topological term: | 2 ,^= 0 ) I = 1- 

This procedure also applies to the SPT internal gauge 
field theory in any other dimensions. We know that 
the SPT state has no intrinsic topological order and the 
SPT’s GSD=1 on any compact closed spatial manifold. 
Therefore, this GSD computation serves as the consis¬ 
tency check that the internal field theory shows a gapped 
phase with nontrivial symmetry transformation — the 
internal gauge field theory realizes SPT state. 

We emphasize that knowing the field theory action is 
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not enough to fully understand the SPT field theory. We 
stress that defining the partition function Z using field 
theory path integral is necessary to fully understand the 
SPT field theory. Below we especially remark the global 


constraints of fields in order to define the SPT path inte¬ 
gral. The partition function in terms of the path integral 
form with a total spacetime dimension d is 


Z = j [Db][Da] exp ( j (^—b I /\da I + 


i(-l ) d - x C 


ijk... , j 


b 1 Ab J Ab K A...)), 


(55) 


here I,J,K , ••• £ {1,2,3,..., d}. Here b is 1-form, a 
is (d — 2)-form, and / = da is (d — l)-form. In the 
presence of symmetry-twist semi-classical background 1- 


form gauge field A, we can write the partition function 
Z as 


Z = J[Db)[Da] exp ( J {-^-{b 1 — A 1 ) A da 1 + ————a b J A b K A ... )) 

= J[Db][Df\exp{J (^( b 1 - A 1 ) A f + l( ~ 1)d ^ ' IJK b 1 A b J A b K A ■..)), (56) 


with the field strength of charges / = da. Importantly, 
we view b 1 and a 1 all dynamical internal gauge fields , so 
they are involved in the path integral measure. 


For f[Df]ef with 


//( 27t) £ Z or Zjv, we 


obtain an analogous constraint on a ID loop: 


Now let us define this path integral properly. Let us 
impose the constraints for this field function in the path 
integral, based on the dual equivalent theory using the 
non-linear er-model. We recall that the a is related to 
the current density j specified by the 17(1) or Zjy charge, 
where we have the total number of charges quantized: 

* J = ffl ^ = ffl ^ *' 27r ' ) € Z ’ *' 57 ^ 


The current density *j is a (d — l)-form, thus jj) of da 

represents the surface integral of a (d — l)-closed mani¬ 
fold, such as a 1-surface for 1+1D spacetime, 2-surface 
for 2+1D spacetime. 

Now we integrate over the field variable / for the parti¬ 
tion function Eq.(56), which procedure analogous to the 
discrete Fourier summation yields a constraint: 

e llpn = 8{ip mod 27r). (58) 

nez 



A 1 ) = 0 mod 27 t 


(59) 


=> j) b 1 = j) A 1 mod 27 t = mod 27r. (60) 

The first line constraint is true for both 17(1) charge 
and Zn charge. The second line constraint Eq.(60) with 

n/ £ Z is an additional constrain if (jfj) f/(2n) 

for our case of discrete Zjy charge for SPT state with 
27v-symmetry. We can still view 6-field sa a 17(1) con¬ 
nection but with a constraint from the Zn symmetry- 
twist probed-field A. This means that the internal gauge 
field b is subject to the global constraint from the semi- 
classical symmetry-twist probed field A. After integrat¬ 
ing out the /, the partition function Eq.(56) subject to 
the global constraint Eq.(60) of the symmetry-twist fields 
A becomes 


[Db] exp (J 

Thus so far by using SPT internal gauge theory path integral, we have recovered the SPT invariant of Ref. [19] claimed 
in the item ( ii ). 



1 Ci 1 i 2 ...i d bh A b i 2 A ... A = exp {J ( l( ^ Cflfa ■ Id A h A A 12 A • • • A A 1 *)). (61) 
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Next, without losing generality, let us take 2 +ID SPT as 
an example, with an explicit Cijk = X) 2, >\ NiN n\-, 3 Pm • 
Let us do the explicit partition function calculation on 
the two topologies, a sphere and a torus respectively, 
by comparing the nontrivial class Z to the trivial class 
ZCpih = 0). For each calculation below we will fix a par¬ 
ticular set of n/ for the global constraint Eq.(60). 

The 1st topology: On a spatial sphere S 2 with a 
time loop S 1 , there is only a non-contractible loop along 
the time direction. So there is only a nonzero nj for the 
global constraints in Eq.(60), and other rij must be zeros. 
We have: 


Z 

Z(pm = 0) 


exp 


1 


= 1. 


(62) 


The 2nd topology: On a spacetime T 3 torus, with¬ 
out losing generality, let us assume, A 1 , A 2 , A 3 along 
x, y, t-directions have nontrivial global constraints with 
some generic n\ , ti 2 and 713 . For example, analogous to 
Sec.VB’s setup, we can assume dx p = dx, dx" = dy, 
dx p = dt. 

Z _exp(ip m ^*i 

Z(;pm = 0) “ I 

(63) 

Since for both on a sphere and on a torus, the ab¬ 
solute value of the above, | Z ( M f I=0 ) I > measures the GSD 
ratio between the nontrivial phases and the trivial insula¬ 
tor. Since the trivial insulator has GSD=1 here, all other 
phases have GSD=1, so the pm 7 ^ 0 phase is a generic 
SPT state. 

We thus confirm that the path integral Eq.(55) with 
dynamical variables describes nontrivial type III SPT 
states in 2+1D. The same procedure can be generalized 
to other dimensions, such as Eq.(20) as SPT states in 
1+1D and Eq.(52) as SPT states in 3+1D. The GSD for 
these theories defined by the partition function is 1. The 
procedure works in more general closed topology, we thus 
show the claim in the item (Hi). 

One further extension of our work is to study the dual¬ 
ity [10] between SPT (which is non-topologically ordered) 
and dynamical topological gauge theory (which is topo¬ 
logically ordered). More precisely, we can start from the 
SPT internal gauge theory path integral of Eq.(56) and 
then dynamically gauge the theory to a dynamical topo¬ 
logical gauge theory equivalent to the Dijkgraaf-Witten 
theory [55]. In Appendix.IX, we will outline such a proce¬ 
dure using field theory path integral, and we will propose 
the continuous dynamical topological gauge theory dual to 
the Dijkgraaf-Witten theory with a discrete gauge group. 


= e 


ipm - 


VIII. EDGE THEORY 

The bulk effective field theory can also describe inter¬ 
esting edge physics. For the 1+1D case, by integrating 


out the Lagrange multiplier fields a 1 in Eq.(20), the cor¬ 
responding edge theory takes a very simple form: 


r° = 

^edge 


'-Cu^d 0 ip J , 


(64) 


with scalar fields ip 1 define the gauge transformation 
ip 1 —7 ip 1 — g 1 to cancel the gauge transformation of 
b 1 —7 b 1 + dg 1 . which is nothing but a quantized topo¬ 
logical term for a quantum mechanical system with de¬ 
generate ground states. Such a Berry phase implies the 
following quantization condition: 


r 1 2i _ 1 _ 27rifVi 2 _ 27 t i 
[<P J Xi ~ piiNiN 2 - pnIV 12 ’ 


(65) 


Here N 12 is defined as the least common multiplier (1cm) 
where N 12 = lcm(A/i, JV 2 ) = N 1 N 2 /N 12 . Due to the 
compactification and the quantization constraint, shown 
in Appendix D 3, the symmetry generators are S v 1 = 

e lNl<p «i 2 and S v 2 = e lN2<p . It is straightforward to 
check that S v i(J dt£° dge )S'" / 1 = (/ dt£° dge )+27r-integer, 

so the partition function Z = f Dip 1 Dip 2 e~ f dt Gdg B j s 
invariant under the symmetry transformation S,„t. We 
find that the symmetry is realized in a projective repre¬ 
sentation manner on the 0D edge, because the symmetry 
generators do not commute: 


_ 2-rvipn 

S v iS v 2 =e w i 2 


( 66 ) 


Here pu is defined as a p\\ (mod N\ 2 ) variable. If 
gcd(pn,Ad 2 ) = 1, it is the Zn 12 Heisenberg algebra and 
requires a Ai 2 -dimensional representation for the sym¬ 
metry generators SLi and S v 2 . This implies the 0+1D 
edge mode of the ground state has a IVi 2 -fbld degeneracy, 
consistent with the edge mode physics analysis via the 
dimensional reduction approach in Ref. [45]. In general, 
even if gcd(p n , N 12 ) + 1, we have a generic gcd( ^]^ 12) - 
dimensional representation for the symmetry generators, 
thus the zero mode degeneracy is 


GSD 


iVl2 

gcd(pn, N 12 )' 


(67) 


For the 2+1D bulk system with its 1+1D edge theory, 
we have an analogous derivation as follows. Integrating 
out a^ leads to the constraint: 

£^ x dX = 0, (68) 


The constraint can be solved by requiring: 


A 


1 

V 


dX, 


(69) 


We see that C c r is nothing but a total derivative: 

£ 0 ff = ^Ci JK £ ll ' 2X d ll ip I d v ip J d x ip K , (70) 
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which actually describes a 1 + ID edge with effective ac¬ 
tion: 

ge = \c IJK e^^d^ J d^ K . (71) 

The higher dimensional generalization is also straight¬ 
forward, e.g., the type-IV SPT in 3+1D can have a 2+1D 
edge theory described by: 

£edge = '-CuKLS^d^d^d^. (72) 

The gapless nature of these boundary terms can be 
proved via dimension reduction to the 1+1D case we dis¬ 
cussed at the beginning of this section. Finally, we note 
that if we view ip 1 as scaling dimension zero fields, £g dge 
and £ 2 dge can be regarded as a fractionalized version of 
0(3) and 0(4) topological theta terms. For future work, 
it would be of great interest to understand the underlying 
conformal field theory described by these fractionalized 
theta terms. 


there exists a duality [ 10] between SPT (which is non- 
topologically ordered) and dynamical topological gauge 
theory (which is topologically ordered). More precisely, 
we can start from the SPT internal gauge theory path 
integral of Eq.(56) and then dynamically gauge the the¬ 
ory by coupling the SPT matter field to external probed 
fields A, and make the A dynamical gauge fields. This 
procedure of gauging SPT with a finite symmetry group 
in principle yields a dynamical topological gauge theory 
equivalent to the Dijkgraaf-Witten theory [55]. Here we 
describe such a procedure using field theory path inte¬ 
gral, and we propose some continuous dynamical topo¬ 
logical gauge theory dual to the Dijkgraaf- Witten theory 
with a discrete gauge group. 


IX. TOPOLOGICAL FIELD THEORY FOR 
DIJKGRAAF-WITTEN LATTICE MODEL 

In Sec. VII, we had established the SPT field theory 
by defining the SPT path integral. It is known that 


Naively, one approach is starting from the path integral 
Eq.(56), if we promote the semi-classical probed field A to 
a dynamical field by including the path integral measure 
[DA], we obtain: 


Z = f [Db][Da][DA] exp ( f (^-(5 J ~ A da J + b 1 A b J A b K A ...)) 


N 


(73) 


One can see that if A is still subject to some global 
constraint: 


, T , „ 27rm , „ 

A mod 27 T =- mod 27 r. 

Ni 


(74) 


but now nj £ Z^ needs not to be fixed. The dynami¬ 
cal gauge theory of A would sum over all possible nj. If 
we compute the GSD of this field theory on a spacetime 
manifold, then we essentially reproduce the same calcula¬ 
tion using the group cohomology cocycle while summing 
over all possible group elements m £ Z/v, . Eq.(73) can 
produce the same physical observables such as GSD of 
Dijkgraaf-Witten theory. This suggests that Eq.(73) can 
be an equivalent description of Dijkgraaf-Witten theory. 

Another approach to obtain the dynamical gauge the¬ 
ory is through the minimal coupling the internal gauge 
field a to the external gauge field A, and then integrating 
out all the internal gauge fields a and b. We describe it 
below. 

2+1D: Now, we are ready to discuss the bulk response 
theory. The external probe gauge field A^ will couple to 


the internal charge current in a standard way: 


Cc 


upling = iApjj = 


(75) 


However, since A 1 is in the Higgs phase with Zm i charge 
condensation, we need to introduce a BF term [42] for 
response guage field A 1 as well: 




(76) 


Actually such a term is crucial for maintaining the gauge 
invariance for the total action.(It is easy to check that 
^coupling is not gauge invariant under the gauge trans¬ 
formation of a 1 ., and we need to shift Bj. to restore the 

r 1 _ h 1 

gauge invariance.) Finally, by integrating out the inter¬ 
nal gauge field a * and , we end up with an effective 
action ig-e^BfaAl + 


_ ^L £ ^x B ig a i , muNiN 2 N 3 x i 2 3 

- 2^ i W 4 «' + (27t) 2 IV 1 93 £ 

( 77 ) 


A- 
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If we view A 1 as background gauge fields describing the 
symmetry twists on the boundary, the above action is 
equivalent to the SPT invariants Eq. (37). However, 
if we view both A 1 and B 1 as dynamical gauge fields, 
the above action potentially describes non-Abelian Berry 
phases, though the original global symmetry is Abelian 
and all the gauge fields are Abelian in its own sec¬ 
tors. The whole Lie algebra becomes non-Abelian fea¬ 
ture due to the central extension Eq.(C4). It will be in¬ 
teresting to verify whether the fully-dynamical topolog¬ 
ical gauge theory is equivalent to the Dijkgraaf-Witten 
gauge theory [55]. Our word of caution is that the non- 
semi-simple Lie algebra detailed in Appendix C suggests 
a more conservative side of this claim. It is also likely 
that method beyond the-saddle-point approximation is 
required to capture the global constraints and missing 
pieces that we may omit in Eqs.(42) and (45). 

3+1D: Similarly, we can discuss the bulk response the¬ 
ory. The external probe gauge field A 1 will couple to the 
internal charge current in a standard way: 

^coupling = i A^fi = -^e^ pa A^a 1 ^, (78) 

Similar to the 2+1D case, we also need to introduce a BF 
term to describe the Z^ I external gauge field in 3+1D: 

^-^B^Al, (79) 

By integrating out the internal gauge field aj tv and A^, 
we end up with an effective action: 

r — tdI a aI 

•‘-response — ^ c ±J p.i/ u p- n a 

-^^4^ + .... (80) 

We warn the reader that there is a potential danger to 
view Eq.(80) as the dynamical topological gauge the¬ 
ory, as one needs to further confirm the physical prop¬ 
erties such as topological GSD and braiding statistics 
must match with the 3+1D Dijkgraaf-Witten topologi¬ 
cal gauge theory [55] computed in Ref. [56]. We will leave 
the study of topological gauge theories for future work. 
The minimum claim of our approach is that viewing the 
B field as a Lagrangian multiplier constrains the flatness 
of A with dA = 0, we essentially derive the SPT invariant 
in terms of the semi-classical probed field A agreed with 
[19]. This confirms our multi-kink topological term and 
vortex condensation mechanism do generate nontrivial 
SPT states. 


X. CONCLUSIONS AND DISCUSSIONS 

In conclusion, we have discussed the multi-kink topo¬ 
logical term and vortex condensation mechanism for 
bosonic Abelian SPT states that cannot be described by 
Abelian Chern-Simons/BF actions. We have pointed out 


that nontrivial SPT states can be viewed as certain Higgs 
phases via defects proliferating in various nontrivial ways. 
Thus, the formalism and concepts developed in this pa¬ 
per can provide further insights for understanding the 
universal mechanism for bosonic SPT states, especially 
for those protected by non-Abelian symmetry. 

Moreover, the general concept of “hydrodynamical ap¬ 
proach” is applicable for fermion systems as well, if the 
spin-manifold is taken into account. Just like we can use 
the spin Chern-Simons theory to describe certain spe¬ 
cial Abelian fermionic SPT states[42], the bulk effective 
actions beyond Chern-Simons/BF theory proposed here 
should also have their corresponding “spin” version that 
can describe new classes of fermion Abelian SPT states. 

The field theory based on the saddle-point approxi¬ 
mation (detailed in Appendix C) may or may not fully 
capture the topological properties of the gapped SPT 
state. However, in Sec. VII, we show that at least for 
the level-1 trivial class of our theory, it has GSD=1 on a 
compact closed manifold just like the SPT state. More¬ 
over, so far as the SPT invariant is concerned, we con¬ 
firm that the bulk SPT response theory induced by the 
multi-kink topological term does reproduce the desired 
SPT invariant. Even though our theory exhibits the so- 
called symmetric-self dual non-semi simple Lie algebra 
[58]; however, due to the extra set of global constraints: 
Eqs.(57), (60), our theory is not equivalent to the usual 
gauge theory with non-semi simple Lie algebra studied 
in the high energy literature (see Appendix C). We be¬ 
lieve our theory is unitary and has finite ground state 
degeneracy on a closed manifold. 

Another important research direction is to study the 
phase transition between superfluids and SPT states, 
analogous to the usual case where we have superfluid 
and insulator phase transition. We will leave these fur¬ 
ther developments for future work. 
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Appendix A: Disorder a superfluid state into a Mott 
insulator or an SPT state 


To guide the readers understanding our formalism, 
here we briefly review this approach using field theory 
(see the pioneer work[47-49] and Ref. [50, 51] for a field 
theory approach). We plan to study SPT states for a 
discrete Abelian symmetry group. First, we will embed 
our discrete Abelian symmetry group into the symmetry 
group of several 17(1) symmetries. Instead of starting 
with a discrete-symmetry breaking state, we will start 
with a symmetry breaking state that break several 17(1) 
symmetries. When we restore the 17(1) symmetries, we 
also restore our real discrete symmetry. 

The superfluid state (the 17(1) symmetry breaking 
state) in any d-spacetime dimension is described by a 
bosonic 17(1) quantum phase kinetic term, whose the par¬ 
tition function Z is: 

Z = J [D9] exp(— J d d x |(0„0 8 + 8 M 0 V ) 2 ) (Al) 

with a smooth piece 9 S and a singular piece 9 V for the 
bosonic phase, and the superfluid compressibility %. We 
stress that the 9 V is essential to capture the vortex core. 
We can introduce an auxiliary field j^ and implement the 
Hubbard-Stratonovich technique [50] , 


Z = J [D9] [Dj^] exp(— J d d z^0?) 2 -iiWs+W)- 

X (A2) 

By integrating out the smooth part J[D6 S ], we obtain a 
constraint 8{d ll j^ 1 ) in the measure of the path integral. 
We can define a generic form 




_ f W! -Vdfl „ 

2n(d-2)\ Wn-w 


with an anti-symmetric a and the total spacetime dimen¬ 
sion d , to satisfy this constraint. More conveniently, in 
the differential form notation, the constraint is d(*j) = 0 
and the resolution is j = 2 ^(*da) with * the Hodge star, 


with an a gauge field in real values. To disorder the su¬ 
perfluid, we have to make the 0-angle strongly fluctuates 
namely we should take the x < Xc or % —> 0 limit 
[51] to achieve large (<9 M 0) 2 . We will however drop the 
Maxwell term due to its irrelevancy in the renormaliza¬ 
tion group (RG) sense. The partition function becomes: 
Z = j\D9 v ][Da] exp(i f d-aA(d 2 9 v )). Hereafter we com¬ 
pensate the dropped ±-sign by redefining the fields. Even 
though naively d 2 = 0, due to the singularity core of 
9 V , the (d 2 0 v ) can be nonzero. Thus, (d 2 0 v ) describes 
the vortex core density and the vortex current, which we 
shall denote (d 2 0 v )/(27r) = *j vortex . In addition, the ac¬ 
tion has a symmetry of a —> a + d£, or more explicitly 
+ % 3 ^ 4 -w]- B 7 Noether theorem, 
this symmetry leads to the conservation of the vortex 
current: the continuity equation d * jvortex = 0, this im¬ 
plies that 


* Jvortex = (d 2 0 v )/(27r) = d6/(27r) 
for some gauge field b. We can thus define the singular 
part of bosonic phase d0 v = 6 as a 1-form gauge field, to 
describe the vortex core, so 

d0 s + d0 v = d0 s + b. (A3) 

The partition function in the disordered state away from 
the superfluid, now becomes that of an insulator state, 
Z = j[Db\ [Da] exp(^ f bA da) with a topological BF ac¬ 
tion. More explicitly, the path integral formalism shows 

Z = J[Db][Da]ex p(i J 2?r ^^ 9 y_ e W2 "' w 6^ 2 a M3 ,.. M J- 

(A4) 

The Hamiltonian of Eq.(A4) is zero, which describes an 
insulator with an energy gap separating the ground state 
from excitations. It has no intrinsic topological order in 
the sense that it has a unique ground state degeneracy 
(GSD, see Ref. [54], this action is a level-1 BF theory with 
GSD=1). This is known as the mechanism of disorder¬ 
ing the charge while condensing the vortices generates a 
trivial insulator: a Mott insulator without SPT order. 


Appendix B: Derivation of the dynamical effective bulk action of SPT 


In the following, we list some details for deriving the internal field theory of SPT in Sec.VB, specifically for type 
III 2+1D SPT with Z^ 1 x Z^ 2 x Z^ 3 symmetry. We note that up to a total derivative, the tri-kink action Eq.(38) 
can be simplified as: 


Ari-kink = + bl) 2 + '-CuKS^-W'Mbib*) - 39^(8^b«) + b^X] + 4 


Maxwell? 


(Bl) 


Again, we can introduce Hubbard-Stratonovich fields jj to decouple the quadratic term as 


Ari-kink = ~ i didX + i Kil + ^CuKsX-WXtbX) ~ W*Md v 0 J X) + bib J X] + Maxwell, (B2) 
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We further introduce Lagrangian multiplier fields and A^ to decouple the — Cijk^^O 1 d tl {d v 6 J b^) term. We 
have: 


Ari-kink = ~ + ibjj^ + W) + bfiX] ~ 


+ iA£(# - CuKe^dvOib*) + C 


Maxwell 


= \{fi) 2 - + # + CuK^bib f) + i bjj? + \c IJK e^%bib f 


+ X(# - C IJK e^dXbl) + ^Lxwell 
1 

2 ( 


— ^O/) 2 — [j/ + £/ + Cijk £^ vX {b J l ,bx - A tb\)\ + ib^jj + ^Cuk^X^X + i + ^Maxwell; 

(B3) 


Integrating out the 6% fields result in a constraint: d^ [jj + + CuK^ vX {bib x — A ib x )\ = 0. From this con¬ 
straint, we can write the conserved jj = j^£^' uX d u a I x — — CijKE flvX (b'lb x — A ib x ). Finally, we obtain: 


Ari-kink = ±e» v %d v a{ - jC IJK e^%b J X + iC IJK e^%\ J X + \ \^Xa[ ~ Cuks^XX ~ *X) 


+ XX + 


i(K ~ K) KZ^ X dXx + C IJK e^ x (bX - *X ) 


2tt 


1 pb 

SJ 1 -^Maxwell 5 


(B4) 


Integrating out the fields, we end up with: 


Ari-kink = X^XdXx - ^CuKE^bWX + i Cj JKE^bi A'&f + £ 


27T 

1 

2 


—e^ A 5,ai - C IJK E^\bX ~ *X) 


i(A£ - tf) - + C IJK E^\bX ~ Aft 


1 2 


27T 




+ £ 


Maxwell 


= + i c ijk e^%\ j X + kK~ K) 2 


2 n 


+ i(Au - b 1 ) 


u fi u iy U X 


/i'V 17 A 


^ ^ x d„a[-c IJK E^ x (b J X - >^X) +£ b 


Maxwell 


Maxwell ’ 


= ±-^%d v a[ + '-CjjKE^XbX iC UK E^%\X+iC IJK e^ x \l\X + ^~ b *X + ^ 

= ^%d v a[ + \C IJK E^\ AjA^Af + (bj, - A^)(6^ - A^)(6f - Af)) + *(&' " ^) 2 F^LweU. (B5) 

r 


Appendix C: Comments on non-semi-simple Lie 
algebra and topological field theory 

In Section V C, we learn that the saddle point approx¬ 
imation leads us to an intrinsic field theory and a bulk 
dynamical theory with non-semi simple Lie algebra. If 
we write the gauge connection in terms of its gauge field 
components and its generators: 

a“T“ = 6^ 7 + a^;. (Cl) 

(dlT\alT 2 ,alT 3 ) = (b^X,, ^X 2 , ^X 3 ), (C2) 
= (alHlalH*,alH* 3 ). (C3) 

Here a = 1,... ,6 and I = 1,..., 3. 

The corresponding generators H 1 and Xj satisfy: 

[Hi ■ H*j} = [Hf , Xj] = 0; [Xj , Xj] = C IJK H* K , (C4) 


where Cijk serves as the structure constant now. The 
full Lie algebra consists of an Abelian Lie algebra X(X) 
with a central extension by another Abelian Lie algebra 
7i*(H*). Here X{X) contains the set of generators Xj 7 
and contains the set of generators Hf. 

For the specific case of level-1 Chern-Simons theory in 
Sec. VII, we are able to show the GSD=1. However, for 
the general level-fc case, the structure of the phase space 
volume is changed (see for example, Appendix D). This 
Appendix is meant to provide some word of caution to 
prevent us from making a stronger claim that the Chern- 
Simons theory with this non-semi simple Lie algebra is 
exactly the dynamical Dijkgraaf-Witten field theory we 
look for, unless we carefully specify the global constraints 
analogous to Sec. VII. 

The particular type of the non-semi simple Lie algebra 
we derived in Eq.(C4) is in the class of symmetric-self 
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dual Lie algebra [58]. Even if the Killing form K ab de¬ 
generates, we can replace the K ab by an invariant nonde¬ 
generate symmetric bilinear form /C G a , if it satisfies the 
criteria below. 

For a Lie algebra given by [T a , T b ] = f a b c T c , the struc¬ 
ture constant f a b satisfies the Jacobi identity: fbc d fad e + 
fca d fbd e + fab d fcd e = 0. The Killing form as a bilinear 
matrix in the adjoin representation can be determined 
from the structure constant. 

Kab = k(T 0) Tb) = —Tr (T a , T b ) = -J2 faa P hp a . (C5) 

a, (3 

The Killing form is called degenerate, if there exists a 
nonzero generator T' such that k(T',T) = 0 for any T. 

In the Euclidean spacetime, we have a Chern-Simons 
theory: 

L=±e*">!CZ a ,(AXx)d v A?(x) 

+ l -fbc a A a ;{x)A b v {x)A c p {x)]. (C6) 

Even if the Killing form is degenerate, as long as this 
{K- G )ij can be found, the (A G )/j can replace the degen¬ 
erate Killing form to make sense of the Chern-Simons 
theory Eq. (C6) with the symmetric-self dual Lie algebra. 

The (/C g )jj is a symmetric non-degenerate invariant 
bilinear form, constrained by: 

fae i (JC G )bi + fab i (JC G )u = 0 (C7) 

The finite and infinitesimal gauge transformations are: 

A„ -»■ A” = U~ 1 {A l _ l + d„)U = + d M )e“° T “, 

A°{x) -A (A“(x) + f bc a A b ll {x)a c {x) + d tl a a {x)). (C8) 

The Lie algebra we find out in Section V C is a sub¬ 
algebra of the most generic symmetric-self dual Lie alge¬ 
bra [58]: 


[X a ,x b ] 

l J ab X + l Jab 

(C9) 

[H a ,H b ] 

— if( h ^ c H 

(CIO) 

[X a ,H b \ 

= if { a f )c X c , 

(Cll) 

[- Ha,m } 

= 

(C12) 


= {H* a ,H*p\ = 0. 

(C13) 


Notice that the subalgebra spanned by X{X) and 
is the Abelian extension of X(X) by 7~L*(H*). 
The full algebra is the semidirect product of 7i(H) by 
this Abelian extension. The particular non-semi-simple 
symmetric-self dual Lie algebra in Eq.(C4) is nilpotent, 
non-abelian, non-reductive and solvable. The corre¬ 
sponding Lie group is non-compact. 

Our theory in Appendix VII is a special case such that 
the GSD is still 1 which can describe the gapped SPT. 
Due to the non-compact Lie group, however, it is likely 


the generic gauge theory of symmetric-self-dual Lie alge¬ 
bra can capture an infinite degenerate gapless phase in¬ 
stead of a phase with finite topological degenerate ground 
states. The concern of (non-)unitarity has been investi¬ 
gated, for example, in Ref.57. 

We believe that the generic difference between our SPT 
path integral and the usual non-semi-simple-Lie-algebra 
gauge theory is the set of global constraints: Eqs.(??), 
(57), (60). For our SPT path integral, the global con¬ 
straints lead to the finite ground state degeneracy, for 
the usual non-semi-simple-Lie-algebra gauge theory, the 
ground state degeneracy can be infinite. It is possi¬ 
ble a more generic theory can describe a state close to 
the potential gapless phase transition between superflu¬ 
ids, symmetry-breaking states and SPT/topologically or¬ 
dered states. We will leave the further investigation open 
for future work. 

Appendix D: Counting the degenerate zero modes 

1. GSD for a gapped system with a 0+1D 
topological term 

We first review a simple ground state degeneracy 
(GSD) calculation by counting the zero mode for a O+ID 
system. Namely we will count the volume of the phase 
space volume, 

GSD = the volume of the phase space, (Dl) 

up to some normalization factor. 

The first system we consider is described by a Berry 
phase term C° = XP. On one hand, in the path integral 
formalism, we have a partition function: 

Z = J [DX][DP] exp[i k J XP], (D2) 

X = OqX is the time derivative X. 

On the other hand, in the quantum operator formal- 

r\ /* 0 

ism, we have the commutator [ X , ^-] = i: 

[X,P] = i (D3) 

X and P are some matrix operators acting on the 0+1D 
space. Here we will consider a compact phase space, 
so that the phase space volume is finite. In particular, 
without losing generality, the identification we assume is 
X ~ X + 27t and P ~ P + 1. Since the Hamiltonian 
is essentially H = XP — C° = 0, the system seems to 
be trivial without kinetic terms or potential terms. How¬ 
ever, there can be degenerated ground states. All ground 
states ’L satisfy H^f = 0. But these ’L may not be all 
independent. To count the GSD thus to count the in¬ 
dependent degree of freedom, we can construct a generic 
ground state VP in terms of the function of A' if we choose 
X as the basis: 

*(X) = J2 c ne inX (D4) 

n£Z 
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The form is obtained by satisfying the constraint: 
V(X) = T(X + 2tt) as X ~ X + 2tt. The 2 tt shift in 
the exponent will not affect the form of the 'k(X) func¬ 
tion. On the other hand, by doing the Fourier transfor¬ 
mation, we can transform the X basis to the P basis via 
'F(P) = f e lfcPA T(A') dP. Up to some normalization 
factor, this yields, 

*GP) = E C n S(kP + n) (D5) 

n£Z 

Meanwhile, the form satisfies the constraint: 'F(P) = 
'F(P-l-l) as P ~ P+l. This implies that c n 5{kP + n) = 
c n -kS(kP + k + (n — k)) . This means that 


c n = c n - k (D6) 

with k G Z. The volume of the phase space is |fc|. 
We have fc independent degenerate ground states de¬ 
termined by k independent coefficients, thus GSD = |fe|. 
The strategy for this example is basically the same as the 
approach in Ref. 54. 


2. Compactification and Quantization 


3. GSD for a gapped system at the O+ID edge of 
1 + 1D SPTs 

After the previous simple first part of calculation, in 
the second part, we consider the 0+1D edge of 1+1D 
SPT. The system we consider is described by a Berry 
phase term in the partition function for the path integral 
formalism: 

Z = j[Dp l }[Dp 2 ] exp[^ J Cu^do^], (D9) 


withC 12 = *^. 

On the other hand, for the canonical quantization with 
quantum opearators, the commutation relation satisfies 




i 27riXi2 

Cl 2 P11N1N2 


(DIO) 


To well-define the denominator for the trivial class pn = 
0, the trivial class’s pu is identified as pu = N\ 2 . 
We may define the conjugate variables as = 

[tp 1 ,C 1 2 '-p 2 ] = i and [p 2 ,P v 2 ] = [ 1 p 2 ,~C 12 p 1 J = i. 

The 1st approach: The compactified size of p 1 and 
p 2 is no larger than 27r, 


For the later convenience, we now set up a relation 
between the constraint of compactification and quantiza¬ 
tion using an angular rotational system as an example, 
with the angle 0 and the angular momentum L. First, 
0 is compactified and identified via: 

0 ~ 0 + 2tt. 

The compactness of 0 leads to the quantization or the 
discretization of its dual variable L, in order to have e iei 
stays invariant as 0 —> 0 + 2-7T. That means, the quanti¬ 
zation is 


A L = 1. 


On the other hand, if we consider the angle 0 is also 
discretized as rotor angle with 


A0 = 


2tt 

iV’ 


then this quantization must come from the compactifica¬ 
tion of L, with 


L ~ L + N. 


p 1 ~ p 1 + 2 tt, p 2 ~ p 2 + 27t. (Dll) 


The quantization and the discreteness of these rotor clock 
is no smaller than: 

* 1 2tt . n 2 tt . , 

A^ = —, Ap = —. (D12) 

Due to the conjugation relation, following the logic of 
Eq.(D7), the compactness in Eq.(Dll) of p 1 ~ p 1 + 2n 
leads to AP^i = C^Ap 2 = 1. Similarly, the compact¬ 
ness of p 2 leads to AP V 2 = C^Ap 1 = 1. Namely, the 
quantization can be: 


27rAi2 ^ 2 __ 2 ttAi2 
Pi\N\ N 2 1 ^ P11N1N2 


(D13) 


On the other hand, following the logic of Eq.(D8), the 
quantization Eq.(D12) implies the possible compactness 
size of P v i and P v 2 as: P v 1 ~ P v 1 + N\ and P v 2 ~ 
P v 2 + N 2 ■ namely, 


<P 


1 


<P 1 + 


2itN 12 
PllN! ’ 




v 2 + 


27 tN 12 

PiiN 2 


(D14) 


In short, due to the constraint of compactification and 
quantization, we have a set of relations: 

0~0 + 2tt^AL = 1, (D7) 

27r 

L~L + No A0=—. (D8) 

The volume of the phase space is A. It can be counted in 
0-space as well as in L-space as (27 t/A0) = (N/AL) = 
N. 


To construct the refined phase space, we need to take 
the largest quantization size in the discretized lattice 
among Eq.(D12) and Eq.(D13), and the smallest com¬ 
pactification size among Eq.(Dll) and Eq.(D14). This 
means that we will require Eq.(D12) and Eq.(D14): 


A^ 1 



27T 

N~2 




1 




2itN 12 

P 11 N 1 




P 2 + 


2ttAi2 

P11N2 
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Therefore the phase-space-volume counting from both 
^j 1 -space and its dual space, </? 2 -space, is both 
^P-/Aw 1 = ip 2 = S™. However, ^ may not 

be integer in general. We will need to multiply a minimal 
factor on the size of the phase space until it becomes an 
integer. This means that in general we will multiply it by 
the minimal phase factor gc<i (p^ jv r ,) until we have an in¬ 
teger size of phase volume: ,, jVl2 . r , = -— ■ v • 

° ^ gcd(pii,ATi 2 ) pn gcd(pn,A'i2) 

The phase-space-volume counting from both (/p^^-space 
and (/? 2 -space results in 


GSD 


Ad 2 

gcd(p n , N 12 )" 


(D15) 


It is straightforward to construct the functional T^ 1 ) 
and its Fourier transformation ^(ip 2 ) with a number of 
gcd(pii 2 ^Vi 2 ) independent coefficients as in Sec.D 1. 

The 2nd approach: We can verify this GSD re¬ 
sult from an alternative viewpoint, by considering the 
projective representation of the symmetry group G = 
Zn x x Zn 2 : 

We propose the symmetry generators as 

S v i =e iNltfi1 ^, (D16) 

=e iN2tp2 ^, (D17) 


state degeneracy at the 0+1D edge, we can study the 
projective representation of the symmetry group acting 
on the zero energy modes, we find: 


S v i S v 2 = e "12 SpiSp 1 . (D18) 

If pu = 0, the symmetry generators are commutative, 
so it can be written as a linear representation; and the 
GSD=1. In general, the symmetry generators are not 
commutative, so it shall be written as a higher dimen¬ 
sional matrix representation. If gcd(pn, Ad 2 ) = 1, it is 
the Z Nl2 Heisenberg algebra and it requires a Ad 2 dimen¬ 
sional representation. This implies the 0+1D edge mode 
of the ground state has GSD = Ad 2 , consistent with the 
edge mode physics analysis via the dimensional reduction 
approach in Ref. [45]. If gcd(pn, Ad 2 ) 7 ^ 1, we can reduce 
the rank of the representation matrix to a smaller rank, 
we rewrite 


iS v 2 = e 


_1 _ pu 

^12 gcd (p II ,M 12 ) 

gcd(p n ,iV 12 ) S V 2S V 1 


In this way, we obtain a relative-prime factor gc< i(p™ n 12 ) > 
and the GSD is the reduced rank of the matrix represen¬ 
tation of the symmetry generators: 


in order to have the symmetry generators invariant under 
the shift over a full compactification size ip 1 —► yd + 27 ^ 3 4 5 

and p 2 — > ip 2 + 2 . Namely, our choice is guaranteed 

to satisfy: S v i(p l ) = S lp i{p 1 + |^) and S^^p 2 ) = 
Sp 2 (p 2 + p^^ 2 2 )• Our choice also obeys the Z Nl and 
Z N2 symmetry: ( S v i) Nl = (S v 2 )^ = 1 when we have 
impose the discretization as Eq.(D12). 

One can check that S<pi(f dt£° dge )S'G 1 = 
(f df£g dge ) + 27 t • integer, so the partition function 

Z = f Dp 1 Dp 2 e~ $ dt£ » d g» is invariant under the sym¬ 
metry transformation S,„r. To calculate the ground 


gcd(pn, ATi 2 ) ’ 

We have shown the degenerate zero modes happening 
on the 0D edge of ID bulk SPT. In general, if we cre¬ 
ate various symmetry breaking domain wall to gap the 
gapless boundary mods of the higher dimensional bound¬ 
aries, we can study the zero modes trapped at the gapped 
domain wall via the dimensional reduction approach. As 
an example, we can look into the 0+1D kink on a 1+1D 
domain-wall edge Eq.(71) of 2+1D bulk SPTs. This re¬ 
sult is consistent with Refs. 19 and 45. 
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